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We give in this letter a set of general rules which allow the prediction of the 

value of the glass transition temperature T g in network glasses. Starting form 

the Gibbs-Di Marzio law which gives a very general relationship between this 

temperature and the average coordination number of a system, we explain how 

to compute from the valencies of the atoms of the glass, the parameter (3 used 

in this law. We check the validity of the obtained expression and show that 

it is possible to predict the glass transition temperature for any composition 

in multicomponent chalcogenide glasses. The possibility of existence of a 

demixed structure is also discussed. 
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Most inorganic solids can be made amorphous by vapor deposition onto cold substrates. 
However, only a very few of inorganic melts can be supercooled by a water or air quench to 
yield bulk glasses which solidify at the glass transition temperature T g . Oxides as vitreous 
silica (S1O2) and chalcogenides (e.g. Ge x Se\- x ) represent some of the best-known glass 
formers in nature. 

The origin of glass-forming tendency and the determination of the value of the glass 
transition temperature is a subject of great interest, not only for purely scientific reasons, 
but also for technological reasons, because the exact or the approximate value of T g is needed 
in some situations in order to optimize the glass preparation process. Therefore, numerous 
efforts have been realized in order to understand the nature of glass transition, but also 
to relate the value of T g to some easily measurable quantities. One of the best known 
relationships is the " two-third'' rule proposed by Kauzmann, stating that T g scales with the 
melting temperature as T g ~ |T m More recently, Tanaka has proposed an empirical 
relationship between T g and the average coordination number. The relationship is readily 
satisfied in oxides, chalcogenides and organic glass-forming materials [Q. 

Importance of thermodynamic factors in glass formation have been discussed by Adam 
and Gibbs ||, and Gibbs and Di Marzio ||, suggesting that the glass transition may be 
a manifestation of a second-order phase transition, when both the Gibbs energy and its 
derivatives remain continous (at the melting temperature, where the liquid crystallize, the 
derivatives of the Gibbs energy are discontinuous, but not the Gibbs energy) . Applying the 
theory to a liquid made of molecular chains, the authors found a quantitative relationship 
between the transition temperature and the density of cross-linking agents inserted inside 
the system. Later on, the relationship was adapted to chalcogenide glasses, yielding the 
following equation, known as the Gibbs-Di Marzio law ||- ||: 

T ° = l-P(<"r>-2) < 01 > 

where T is the glass transition temperature of the chain-like glass (e.g. vitreous selenium 
with To = 316 K), (3 a system-dependent parameter and < r > the average coordination 
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number. The latter has been introduced by Phillips in his constraint theory and is widely 
used in the investigation of network glasses ||- ||10|| . < r > is defined in a M-component 
glass by < r >= Y^=i m i x ii where mi is the valence of an atom with concentration Xi (e.g. 
< r >= 2.67 in GeSe2)- Agreement of the Gibbs-Di Marzio law with experimental data 
could be obtained by fitting the parameter (3 for numerous glass systems |J- |J. Thus, the 
Gibbs-Di Marzio law seems to describe very well the T g trends as a function of the average 
coordination number < r >, at least for concentrations corresponding to < r >< 2.7. 

On the contrary, if one could compute exactly the parameter (3 of a given glass system, 
one could be able to predict its glass transition temperature T g as a function of < r >. The 
purpose of the present letter is to show that the parameter (3 is related to the local glass 
structure and can be easily computed for any chalcogenide glass system. 



In a recent series of papers [|TT[|- [14]], R. Kerner and one of us have demonstrated by use 
of a model of statistical agglomeration that the glass transition temperature T g displayed a 
very simple law (called " slope equations 11 in weakly modified binary glass systems, such as 
Ge x Sei_ x with x the concentration of the modifier atom less than 0.1. 

T 



dT 9 

dx 



x=0,T g =T J n 



(0.2) 
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where To stands for the initial glass transition temperature and has the same meaning as the 
one appearing in equation ( p.l|) . is the coordination number (or valence) of the modifier 
atom (e.g. m# = 4 is the valence of the germanium atom in Ge x Se\_ x systems). The value 
of these coordination numbers can be determined in most of the situations by the 8 — N rule, 
where N is the number of outer shell electrons of the considered atom | 15| . Thus, the slope 



in equation ( p.2|) indicates how the value of the glass transition temperature will change if 
a small proportion of modifier atoms (such as Ge or As) is added to the initial network. 
The relationship ( p.2| ) is very well satisfied in more than 30 different glass systems such as 
chalcogenides and binary glasses (e.g. (1 — x)SiC>2 — xLi 2 0) Last but not least, if one 
uses the average coordination number of the network < r >= mgi + 2(1 — x) and performs 
the first-order Taylor expansion of the Gibbs-Di Marzio law (|0.1| ) in the vicinity of < r >= 2, 



it is possible to obtain an analytical expression for the parameter (3 in two- component glass: 



0: 



1 



{tub — 2) In 



m B 



(0.3) 



Here again, the predicted value of the parameter (3 computed from obvious structural con- 
siderations (rriB = 4, tub = 3, etc.) is in excellent agreement with the value determined 



from experimental data |12| . We shall now prove that the factor appearing in the right-hand 
side of equation (p.3|) has a universal character and can be easily extended to a multicom- 
ponent chalcogenide glass system, yielding the value of the parameter f3 for any sytem, to 
be inserted in the Gibbs-Di Marzio law ( |0.1| ). To do this, we shall compute from available 
experimental data the value of the parameter (3 and compare it to the predicted one. 

Let us extend the expression (p.3|) to a M-component glass. For the reader's convenience, 
we shall first consider a glass system made of three different kinds of atoms (say A, B and 
C, with respective concentration 1 — x — y, x and y), one of them being the atom of the 
chain-like initial structure (when x = and y = 0). The valences of the involved atoms are 
rriA = 2, itib and mc- The average coordination number is < r >= mBX+rricy + 2(1— x—y). 
We can compute the derivative with respect to the glass transition temperature: 



d < r > , . dx^ , . dy 

-ST = ( " lB " ^ + (mc " 



(0.4) 



and look at the limit when < r >~ 2 (i.e. T g ~ To, x = and y — 0). Then, we can still 
identify the derivative of the first-order Taylor expansion of the Gibbs-Di Marzio law (p.l|) 
with the right-hand side of equation (p.4|), where the quantities dx/dT g and dy/dT g have 
the form presented in equation (p.2|) . Identifying all this and simplifying by To leads to the 
analytical expression of the parameter (3 in a glass made of three components: 



- = (m B - 2) In 



m B 



+ (m c - 2) In 



m c 



(0.5) 



The extension to multicomponent systems appears to be quite natural. f3 has the same 
sum rules as the resistance in a parallel circuit in electrokinetics, i.e. it is the sum of the 



1/(3 of each related two- component system AB, AC, etc. For a given system made of M 
different kinds of atoms with valencies m^, we just have to sum up the M — 1 contributions 
(rrii — 2) ln[^.] in order to obtain the theoretical value of f3~ l . 



i = £K-2)ln 



(0.6) 



i=l 

In this notation, tum is the coordination number of the chain atom, equal to 2. We have 
looked at the validity of this expression on a variety of different glass systems (ternary glasses 
as As x Ge y Sei- x - y and multicomponent glasses such as As — Si — Ge — Se, etc.) 

Comparison with experimental data: In order to minimize the influence of the preparation 
techniques and to obtain a meaningful correlation coefficient, we have carefully selected data 
of i) glass systems prepared with the same heating/ cooling rate ii) glass systems with more 
than five different compositions. We have checked that the influence of the heating rate on 
the value of j3 could be neglected at the heating rates which were used in the preparation of 



the glasses we have investigated. To do this, we have used both the Kissinger's formula [16 

In 



T 2 
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+ const. = — (0.7) 



Q 

which is very well adapted for the description of chalcogenide glasses, and the Gibbs-Di 
Marzio law. Q is the heating rate, E the activation energy for glass transition and R the gas 
constant. Inserting typical values of < r >, E and T shows that /3 depends weakly on Q 
(A/3/ (3 < 5%). This verification has been also realized numerically on systems for which the 
heating rate in DSC calorimetry was reported. For example, the parameter (3 of the glass 
As — Sb — Se lies in the range [1.14, 1.17] (computed from a least-squares fit) for cooling 
rates between Q = 0.62 K.min~ l and 40 K.min -1 , which corresponds to A(3//3 = 0.013. 

The initial value T has been averaged over a set of data found in the literature (T 
of v — Se has been taken as 316 K, of v — S as 245 K). We have performed a least- 
squares fit of the Gibbs-Di Marzio law applied to the data. The results of the fit for the 
parameter (3, denoted as /3 exp , and the correlation coefficient are displayed in Table I with 
their corresponding reference. For completeness, we have reported in the first part of Table 



I some results [[12] for two-component systems, which satisfy (|0.3| ). The results can be 



compared with the predicted value of /3, denoted as (3 pr in Table I, which has been computed 
from equation ( p.6| ). For example, in the system Si-As-Ge-Te the involved valencies of the 
modifier atoms are 4, 3 and 4 respectively Thus, /3Z. 1 = 3 In 2 + In 3 and (3 pr = 0.31, in 
excellent agreement with the fit f3 exp = 0.30. The number of different glasses satisfying 
exactly or roughly the rule ( p.6|) is impressive and proves that the agreement is not a matter 
of coincidence. 

Let us now come to the prediction of glass transition temperature values. We can note 
that the addition of a two-valenced atom (as tellurium or sulphur) in multicomponent glass 
systems does not affect the value of the parameter (3. A 4-component system which involves 
a two-valenced atoms can therefore be considered as a ternary system (Table I) and T is 
then the glass transition temperature of the initial mixture (e.g. Se — Te). But for all the 
other elements of the columns III, IV and V, we can compute the value of f3 and represent T g 
as a function of < r > (fig.l). The simultaneous use of equations ( p.l[ ) and ( |0.6| ) should give 
the value of the glass transition temperature of any composition, at least for < r >< 2.4 (as 
shown on figure 1 for the comparison between the theoretical Gibbs-Di Marzio law and the 
experimental data). For greater values of < r >, one should take into account intermediate 
range order effects such as the existence of rings |J- [|l(|. This extra influence has not been 
considered here. 

Influence of a demixed structure: Nevertheless, there are still several systems for which 
the predicted value of (3 does not match with the computed one. The difference between (3 exp 
and (3 pr can be quantitatively discussed in terms of the presence of some demixed structure 
inside the network. This feature is generally used in order to explain the unusual variation 
of T g in glasses, such as Bi x Se\_ x or ln x Se\_ x . In the former system, data obtained for 
the Bi — Se glass, show that glass transition temperature increases with Bi content for low 
modification, in perfect agreement with the slope equation ( p.2| ). But when x > 0.1, T g 



remains constant P21. This behavior is attributed to the existence of a demixed structure 



of Bi2Se^ microclusters. The same happens in the indium selenide compound at every 



concentration, and evidence of In 2 Se^ clusters has also been discussed |23 . 

Again, let us consider as before a ternary glass system A\_ x _ y B x C y with the correspond- 
ing coordination numbers rriA = 2, mg and mc- If we assume that the general tendency of 
the glass is to form a demixed strcuture of A and B in stoichiometric proportions, then we 

can rewrite the system as: (B rnA A mB ) x / mA CyAi_ y _( mA+mB } x / rnA which defines an effective 

y 

l-i(raj+ra B )/ra,4 

the system is given by: 



concentration of C atoms v e ^ = - — -, r — 77 — and the average coordination number of 



(m c -m A )y 

< r > = m A H 7 — t — r (0.8) 

j _ {m A +m B ) 



m A 

If we proceed as before, i.e. taking the derivative with respect to T g and looking at the limit 
(x, y — > 0) in order to identify with the expansion of the Gibbs-Di Marzio law, we can see 
that the corresponding parameter (3 is defined as (3~ l = (mc — m^) ln[^]. In other words, 
the parameter (3 of a ternary system which displays a demixed structure can be computed 
by considering only the remaining two-component glass. Here again, we have checked the 
validity of this rule on a set of germanium incorporated chalcogenides (Table II). We have 
considered all possible stoichiometric demixed structures at the tie-line composition (e.g. in 
Ge x As y Sei- x -y, the possible structures are GeSei and As^Se^). Again, we have computed 
the corresponding parameter {3 pr and compared it to the value f3 exp obtained from a least- 
squares fit. Table II shows that most of the III — IV — VI systems such as Sb — Ge — S 
glasses behave as a single IV — VI glass (as Ge — Se) with parameter j3 pr = 0.72. This is 
explained by the presence of A 2 X 3 (X = S, Se, Te) clusters with A an element of the column 
III or V. 

In summary, we would like to recall the most important rules to be used if one wishes to 
determine a glass transition temperature of a multicomponent chalcogenide glass. 

1 . Consider all the possible demixed stoichiometric structures of the M-component chalco- 
genide glass. There is in general a strong evidence of demixing in the case of Sb ternary 
based glasses 
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2. Given the concentrations Xi of the atoms, compute the average coordination number 
< r >, either directly, or by use of the effective concentration if there is some 
evidence of structural demixing. 

3. Compute the parameter (3 from the coordination numbers of the remaining atoms 
(those which are not involved in the demixed structure). 

4. Insert (3 and < r > in the Gibbs-Di Marzio law in order to obtain the glass transition 
temperature. 

As illustrative examples of these rules, we give for the conclusion glass transition tempera- 
tures of systems, which have to our knowledge never been investigated. 

Assuming that A/3//3 ~ 0.1 (possibly produced by the heating rate effects, see above, or 
anything else) and if there is no demixing, the glass As^Ga^Sii^Se^ should have a transition 
temperature of about 385 ±8 K since (3 = 0.45 is computed from the valencies m 8 = (3, 3, 4) 
and < r >= 2.4. In the sulphide analog glass, one should measure 299 ± 7 K. Similarly, 
the glass Gaio-Bio^so should have a transition temperature of T g = 325 ± 11 K (/3 = 1.23 
computed from the valencies rrii = (3, 3) and < r >= 2.2). 

The influence of intermediate range order on this set of rules will be discussed in a 
forthcoming article. Also, work devoted to the influence of the cooling/heating rate on the 
value of the glass transition temperature is in progress. 

The authors gratefully acknowledge R. Kerner, R.A. Barrio and J. Ledru for stimulating 
discussions. This work has been supported by the CONACyT grant No. 25237-E. 
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System 



Ge-(Se) 
Ge-(S) 
Si-(Se) 



Ge-Sb-(Se) 
Ga-Ge-(Se) 
Ga-Ge-(S) 
As-Sb-(Se) 
Al-P-(Se) 
Ge-Sb-Te-(Se)* 
Si-As-Ge-(Te) 
Ge-Sb-As-Te-(Se) 



P P r Pexp Correlation Reference 
coefficient 



0.72 
0.72 
0.72 



0.72 
0.73 
0.81 



0.988 
0.998 
0.997 



from 12 



from JT2[ 
from [12] 



0.56 
0.56 
0.56 
1.23 
0.32 
0.56 
0.31 
0.45 



0.66 
0.55 
0.59 
1.17 
0.21 
0.55 
0.30 
0.55 



0.972 
0.965 
0.823 
0.995 
0.952 
0.998 
0.979 
0.989 



I 



24 

m 

191 



20 
§ 

§ 



TABLE I. Different multicomponent glass systems. Comparison between the predicted value (5 pr 
obtained from equation ( |0.6| ) given the valencies of the involved atoms, and the value (5 exp computed 
from experimental data by a least-squares fit. The atom which corresponds to the chain-like entity 
corresponds to the chemical symbol inside the bracketts. M-component glasses with an * can be 
considered as (M-l) component glasses for the computation of f3, since they involve a two-valenced 
atom (Te). 
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System 



In-Ge-(Se) 

Sb-Ge-(S) 

Sb-Ge-(Se) 

Sb-Ge-(Te) 

Sn-Ge-(Se) 



(3 exp Demixed Correlation Reference 
structure coefficient 



0.77 In 2 Se 3 

0.61 Sb 2 S 3 

0.78 Sb 2 Se 3 

0.79 Sb 2 Te 3 

0.68 SnSe 2 



0.993 
0.932 
0.986 
0.992 
0.973 



m 
m 
m 

WL 
H 



TABLE II. Different multicomponent glass systems exhibiting a demixed structure. Comparison 
between the predicted value of (3 pr = 0.72 obtained by considering the remaining Ge — X system 
(X = S,Se, Te) and the value (3 exp computed from experimental data by a least-squares fit. The 
atom which corresponds to the chain-like entity corresponds to the chemical symbol inside the 
bracketts. 
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FIG. 1. Several Gibbs-Di Marzio laws with the computed value of (3 in multicomponent chalco- 
genide glasses. Solid line: mi = (3,4,4) and (3 = 0.31. Dashed line: mi = (2,3,4) and (3 = 0.56 
with the experimental data of the Te — Sb — Ge — Se system [3] (boxes). Dotted line: rrii = (3, 3, 3) 
and (3 = 0.82. 
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